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Abstract 

o, 

'nT , Applying the bosonization procedure to constrained fermions in the frame- 
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work of the one dimensional t — J model we discuss a scenario of singlet 
superconductivity in a lightly doped double chain where all spin excitations 



G , remain gapful. 



X 



Last years the problem of non-Landau Fermi liquid behavior in quasi-one dimensional 
systems again attracted a strong interest. This time it was stimulated by the Anderson's 
idea about an effectively one dimensional (ID) dynamics of excitations in the normal state 
of layered high Tc cuprates |I[] . 

The basic problem here is a complete description of a " dimensional crossover" which may 
occur as a result of varying coupling between ID systems (chains) forming a two- or three- 
dimensional array. In particular, one of the main issues is whether a coherent transport 
between chains establishes at arbitrary small interchain coupling or whether there is a finite 
threshold resulting from the "confinement" phenomenon |jl]]. 

Various weak coupling studies of the infinite array problem don't seem to confirm the 
Anderson's picture, although one might think that the situation becomes different at strong 
coupling. On the other hand systems of a finite number of chains provide interesting exam- 
ples of a peculiar behavior in an "intermediate dimenion". In analogy with a purely ID case 
one might expect that these also allow a consistent strong coupling treatment. 

Moreover these models can also describe properties of such real materials as (^0)2^207 
or Sr2Cu/iOQ which contain weakly coupled metal-oxide-metal double chain ladders. It 
was also pointed out in [Q that higher stoichiometric compounds Srn-iCun+i02n provide a 
physical realization of weakly coupled A^-chain ladders. 

A proper Hamiltonian of a strongly correlated double chain is that of the t — J model 
with antiferromanetic spin exchanges 0: 

H = -Y.i^H cia,fCi+l,aJ + t±cl^^^Ci,a,d + h.C.) + 

+ ^{JSijSi+ij + J±Si^uSi4), (0.1) 

i 

where f = u,d is the chain index. The Hamiltonian (1) has to be complemented by the no 
double occupancy constraint 

E 4.jc,.,/ < 1 (0.2) 

Accordingly, a weak coupling regime can be studied in the framework of the Hubbard model. 



At half filling the Hamiltonian (1) describes 5* = 1/2 double chain Heisenberg model. 
Available numerical results 0, Q as well as a mean field analysis using the Gutzwiller 
projection ||^ indicate that in contrast to the case of a single chain there are no gapless spin 
excitations in a double chain. Moreover the spin gap appears to be robust against doping 
and survives in some range around half filling. These observations are in agreement with the 
conjecture that a lightly doped double chain system becomes a singlet superconductor of 
a modified d-wave type. 

Recent weak coupling renormalization group (RG) studies of the double-chain Hubbard 
model did reveal some spin gapful fixed points characterized by an inhanced singlet pairing 
in both cases f/ < t ~ t_|_ [^] and t± < U < t ^. An inhancement and a power-law decay of 
pairing correlations were also convincingly shown numerically Q]. 

To clarify the essence of the double chain physics it is worthwhile to review properties 
of the single chain t — J model. 

In the region of ratios J/t < 1 the model can be only found in the so-called Tomonaga- 
Luttinger (TL) regime which corresponds to both gapless spin and charge excitations [Q, 
0. It is customary to describe the TL behavior in terms of spin and charge correlation 
exponents Kg and Kc- 

The spin exponent Kg equals to unity everywhere in the TL regime while Kc gradually 
increases from the value 1/2 which it reaches at J = and any density p as well as at p ^ 1 
and arbirary J/t as J/t increases or p gets smaller. The TL regime persists up to J/t ^ 2.5 
where a spin gap with strong pairing correlations occurs at small enough fermion density 
(p < 1/3). In fact, one can understand the occurence of the region of attraction at small p as 
resulting from the existence of a two-particle bound state at zero density. A finite threshold 
in the attraction strength follows from vanishing of the bound state wave function at zero 
separation due to the no-double occupancy constraint. 

On the other hand in the regime of strong correlation at p close to unity one can argue 
that gapless spin fluctuations drive couplings of the charge sector to the repulsive region and 
to get an effective attraction {Kc > 1) one has to exceed some threshold value of J/t. 



However if both p and the critical value of J /t are large then the attraction of charges 
actually leads to the phase separation rather than to the real superconducting pairing. Even 
the inclusion of the short-range repulsion which is supposed to postpone the onset of phase 



separation to higher J /t doesn't extend the region of singlet pairing |T0 . 

If, on the contrary, spin fluctuations are gapful then these may not renormalize charge 
couplings significantly. Therefore the charge correlation exponent may not receive its basic 
contribution leading to i^c < 1 at small J /t and the mechanism of attraction may work 
without any threshold in J/t. Various possibilities to get such a behavior by means of the 



frustrating spin exchange interactions in a single chain were considered in [|ri|], ||T2 |. 

Coupling between Luttinger chains provides an alternative way to produce a spin gap 
favoring attraction between charges. 

Preceding analytical studies of the problem in the framework of the bosonized "g-ology" 



0, [|T3|, |0| were restricted on the case of weak coupling. 

In the present paper we find additional arguments in favor of the above scenario by 
using a bosonic representation of the t — J model which is an adequate tool to study a 
strong coupling behavior. 

Although the method of bosonization is conventionally applied to ID weakly interacting 
fermions it might be possible to formulate a consistent procedure for the opposite limit when 
the interaction is extremely strong, that is for the case of constrained fermions. Various 
versions of the bosonization procedure in the framework of the t — J model were discussed 
in the literature |T^, [|1^. Recently a modification of the approach proposed in [|15 



was 



shown to give correct exponents for the one- and two-particle correlation functions as well 
a good approximation for the energy spectrum of the single chain t — J model |jT6[ . 



According to the method of |T6[ the constrained fermion operator Cj o- can be represented 
as a product of a spinless fermion \I^j and a spin one-half operator (hard-core boson) S^ : 

where the projection operator Pi reduces the space of four on-site states {\hole > ®\spin > = 



|1, t>, |1, |>, |0, |>, |0, |>) to the physical Hilbert space formed by the set |0 >, | t>7 1 i>- 
In turn, the spin one-half operator Sf" can be expressed in terms of a spinless Jordan- 
Wigner fermion 

St = xle^^^<^'^''' 5- = x.e-^'^S.-^^^^ S^ = xlx^-l (0.4) 

The authors of [jT6[ also argued that the local constraint (2) and the sum rule for the 
constrained fermions J2a ItT ^^"^ ^ {4o-('^)5 Ci,<j(— ti;)} >= 1 + <5 (where 5 is doping) are 
obeyed even in the approximation which discards the projector Pi. Therefore this approx- 
imation which reportedly yields correct exponents of correlation functions and momentum 
distribution was supposed to provide a better account of fermion correlations than other 
approaches which don't treat hard-core condition properly. 

With the neglect of projectors Pi the representation (3) allows one to rewrite (1) in the 
form 

H = - E(iE^l/^«+ij(^^/^r+i,/ + SrfSt,,j) + tM,J!^4iStuS■,, + Sr^S^,) + h.c.) + 

i,u f 

+ Y^iJ^ijSijSi+ijTii+ij + J±ni^uSi,uSi,dni^d) - fJ'Yl ^l,/^« J (0-5) 

with a local charge density defined as riij = \E^| yVE'jj and a local spin Sij given in terms of 
Xij according to (4). 

To get a bosonic form of the Hamiltonian (5) one can use the representation 

fi=R,L fi=R,L 

Keeping the most relevant operators in the continuous limit we obtain the bosonized Hamil- 
tonian density 

+vs{{deff + {dct^^f) + J{l-5 + -^{dct^t ± d^c)? cos2(C ± 0+)} + 

-|-t_L(cos0~ + cos(0^ + 25x) cos0^) cos6'^ cos6'g' + 

+Jx(l - 5 - -^{dcPt + d<p-)){l - 5 - -^{d<Pt - aC))(cos20+ + cos 20; + cos^;) (0.7) 

v27r yzTT 
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where 0^^ = ^^(^"s ^ 4>i,s) ^^^ Of,, = -j^{Oc,s ^ ^c,s)- -^^ usual, the values of correlation 
exponents Kf^ can be affected by short-wavelength renormalizations. The bare velocities of 
charge and spin excitations are given by the formulae: Vc = 2t sin nS and Vg = 2 J((l — 6)"^ — 
(^^^)^) + 47rtsin7r5 although these can be altered too. 

We note that the expression (7) can be also obtained by using the CP^ coherent state 
representation first introduced in ll^ and applying a somewhat different bosonization scheme 



discussed in [|^. We believe that it verifies a neglect of the projecting operators when 
deriving a relevant part of the continuous bosonic Hamiltonian. 

The first three terms in (7) can be recognised as the Hamiltonian of the charge carry- 
ing spinless fermion coupled to the Abelian gauge field A^ = e^j^udyCps = {d(l)s,d9s) which 
describes a surrounding spin background. Notice that once the constraint was explicitly 
resolved one obtains only two independent fields instead of three (spinon, holon and a gauge 
field) appearing in mean field studies of the t — J model. 

In the case of a single chain the basic TL phase of the t — J model corresponds to both 
(pc and (ps being gapless. Scaling dimensions of operators of the form cos /50 cos j3'6 can be 
estimated by means of the formula 

A = \{K(3' + ^) (0.8) 

At (5 = (7) becomes equivalent to the bosonized Heisenberg double chain due to the 
effective freezing of the charge degrees of freedom. The resulting expression essentially 



coincides with the relevant part of the one obtained in ||T8| where a more general XXZ 
symmetrical case was considered. 



It was arued in |]I8[ that at least one of the operators cos 20~ and cos 9j appears to be 
relevant and drives the "-" spin sector toward a strong coupling regime where either (j)~ or 
6j gets locked and a corresponding cosine acquires a nonzero expectation value. 

Additionally, at Kf < 1 the "+" sector gets to a strong coupling regime where (pf 
is locked. Therefore in general there are four possible phases with finite Vs which can be 
identified as follows {rj = 2/fC+): 
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-Both 0~, 0+ are locked- the antiferromagnetically ordered state ; 

- Both 6~ and cj)^ are locked - singlet state (all spin excitations are gapful); 
-9~ is locked - the XY-type phase characterised by correlations: 

<S}{x)S}{^)>r^^ + {-ire-^ 
< Slix)SjiO) >~ (-1)-^ + e-- (0.9) 

- (f)~ is locked - another gapless phase having different spin correlations: 



<5/(x)S7(0) >~e-^ (0.10) 



All these states besides the last one were argued in [|T8[ to appear on an extended phase 



diagram of the double chain XX Z symmetrical (Heisenberg-Ising) spin model. The analysis 
carried out in [|18| leads to the conclusion that at J_l > — | J^ the gapless line J± = becomes 



unstable against arbitrary small Jj_ of any sign. 

It has to be noticed that at J_l < and J± > the nature of the singlet ground state 
is quite different. In the former case every pair of spins on one rung of the ladder tends to 
form a S = 1 state and the system effectively behaves similar to S* = 1 Heisenberg chain 
while in the latter case spin pairs couple preferably into singlets which then form a "dimer 
liquid" . 

According to the phase diagram proposed in [1^ the spin isotropic poin at J^ > is 



located deeply inside the gapful " dimer liquid" phase with both 9^ , 0+ being locked. It 
agrees with numerical 0, 0] and mean field results. 

We also note that in a general case of A^-chain ladders one may expect that a spin gap is 
present at even A^ only. To see that one can apply arguments due to Haldane [|19|. In |]19| a 
topological term governing a longwavelength dynamics of the 2D lattice Heisenberg model 
was found in the form J2y{—^yQy{x, t), where Qy{x, t) is a topological 6'-term appearing in a 
purely ID case and distinguishing between integer and half-odd integer spins [^. In the 2D 



case with periodic boundary conditions this sum is equal to zero which means the absence 



of a 2D counterpart of the ID 0-term. However applying this formula to the finite width 
strip one can see that for odd A^ the above sum doesn't vanish and therefore an effectively 
ID longwavelength dynamics remains gapless. In contrast, for even N the gap survives and 
then scales as A(A^) ~ exp(— A^). 

On the basis of the results obtained in |]^ , ||^ , [^] we assume that at small doping the 
hopping terms in (7) can be treated as perturbations which do not destroy the gapful spin 
state. Then we get the effective Hamiltonian describing a charge dynamics in the spin gap 
state 

Hc = Iy: vf{Kf{d9ff + -^(90^)2) + h COS0,- cos^: (0.11) 

where i± = t_j_ < cos 6 J >. Charge correlation exponents which can be easily read off from 
(5) are given by the formulae 

Kf = {! + -< SfSf > ±4 < SuSd >)-'/' (0.12) 

TTt ■' ^ Tit 

Due to the short-range antiferromagnetic order we encounter the case of K^ > 1 while 
K^ — 1 can be, in principal, of both signs. 

The physical origin of attraction between charges in the paramagnetic spin gap state 



with a short range antiferromagnetic order can be understood on very general grounds ^1 
A straightforward manifestation of this phenomenon in the framework of the t — J model is a 
negative sign in front of the product of charge densities nfUfi staying in (5) if < SfSf >< 0. 



A progressive understanding of the spinless double chain problem |22|, ||23|, 1^^ shows 
that despite of the possibly vanishing of the single particle hopping the "-" charge sector of 
the system always evolves to the strong coupling regime due to the development of either 
coherent particle-particle or particle-hole pair hopping. This phenomenon was previously 
discussed by many authors in the context of quasi-one dimensional conductors pSf . 

For the account of these processes triggered by the single particle hopping the charge 
Hamiltonian (11) has to be supplemented by the extra terms 

SHc = Qph cos 20~ + Qpp cos 26~ (0. 13) 



generated in the course of renormalization. At small t^/t a conventional RG procedure 
applied to the extended Hamiltonian (11,13) leads to the system of equations describing a 
renormalization flow in the "-" charge sector |2^, [J^ (■C is a scaling variable): 



dgph 

d^ 

dQpp 

~dl 



2{l-K;)g,^ + ti{K 



K: 



-) 



2{l-^)gpp + ii{^-K;) 

d\ogK~ 1 2 , 1 2^ 



(0.14) 



di 2' 

The analysis of the solutions of (14) first performed in p4| shows that depending on the 
sign of g- = gph — gpp either cos 20^ (at (7„ > 0) or cos 29^ (at g- < 0) acquires a nonzero 
expectation value. The asymptotic behaviors of the correlation exponent in the two cases 
are K~{^) -^ and K^{C,) — > oo respectively. By considering four possible order parameters 
for the spinless case 



CDW+ ~ cos( 



SS+ ~ cos(0+ + 0-^ 



CDW- ~ cos( 



SS^ ~ cos(^+ + 



e: 



C ' TC 



(0.15) 



one concludes that at (?_ > the competing types of ordering are (intrachain) CDW^ and 
(interchain) SS^ while a.t g^ < the relevant orderings are CDW- and SSj^. In turn, the 
result of the competition between them depends on the sign of (7+ = gph + gpp- 

By mapping the Hamiltonian (11,13) onto the spin S = 1/2 chain in an external magnetic 



field the authors of |2J] also argued that the above statements hold at strong coupling too. 
Considering different order parameters relevant for fermions with spin 






COS 



cos(0+ + 



CDW+ = Y, 

CDW. = y: ciiczi^ ~ cos(0+ + e:) cos(0+ + e. 



SDW+ = Y 



J] J 



c,'c 



cosidij + 



cos 



SDW. = Y: iicZi^^ ~ cos(0+ + d:) cos(e + C 



ss^ = Yj ^c 



iaC-f^,-a ~ C0S( 



^c + fc ) ^m<Ps 



TS+ = Y: ^clci,,. ~ cos(e + 0-) sm(e + o;) 

TS^ = Y: crci^cL,,^ ~ cos(e + 0,") sm(e + 07), (0.16) 

we observe that in contrast to the spinless case, any paramagnetic spin gapful state has only 
two relevant types of orderings which are SS- and CDW- containing fields 6~ and (pf ^. 

Examining the divergencies of corresponding response functions with the use of (8) and 
(12) we obtain that at J± < SuSd >< the interchain singlet pairing SS^ appears to be 
the leading instability while in the opposite case J± < S^Sd » the ground state is the 
"fiux phase" CDW_. 

The latter state is characterised by the commensurate with density "fiux" $ = 2kF 
defined as a circulation of a phase of the on-rung order parameter < uldi + d\ui > through 
a plaquette formed by two adjacent rungs of the ladder. In the case of spinless fermions 



cLS cL 



this type of ordering called the " Orbital Antiferromagnet" was first discovered in |2^ 
counterpart of 2D fiux states. 

It is also instructive to express the above order parameters in terms of the hybridized 
states corresponding to the mean field "bonding" and " antibonding" bands B,A= -j^{u±d): 

ss^ = Y 4.4,-. - <4,-. (0.17) 

Considering the distribution of signs of the order parameter SS- on the "four-point Fermi 
surface" {k = (kpyO), {—kF,0), {kp,7i), {—kp,7i)) we observe that it corresponds to the "d- 
wave" type pairing. One might expect that in a two-dimensional array of weakly coupled 
double chains with a continuous Fermi surface this type of ordering does transform into an 
ordinary d-wave pairing. 

It follows from the preceding discussion that both instabilities develop without a thresh- 
old in J/t or J±/t. Note that this statement is in agreement with the results of the weak 
coupling analysis of the double chain Hubbard model . 
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We also add that in the antiferroinagnetically ordered state where all spin excitations 
acquire Ising gaps, the only relevant order parameters could be SS^ and CDW^ containing 
fields (f)~ and 0^. However at all J < SfSj > < Jj_/t the intrachain pairing SS+ is always 
favored. 

In summary, in the present paper we applied the method of bosonization to constrained 
fermions in the context of the double chain t — J model. As a result, we found further 
arguments supporting the recently proposed scenario of singlet superconductivity in the 
spin gap state of the double chain problem |^. 

Our analysis was based on the assumption that a small doping doesn't destroy the 
spin gap and a spin dynamics remains essentially the same as in the insulating case. This 
conjecture is strongly supported by the results of numerical studies [§], 0], the Gutzwiller 
projected mean field 0] and weak coupling "g-ology" P, 0|- 

We also want to stress that our conclusions contradict with a recent claim p^ about an 



existence of the strong coupling fixed point where some spin excitations remain gapless. The 
authors of P^ considered the double chain t — J model without an interchain spin exchange 



(J^ = 0). Then on the bare level their Hamiltonian can be assigned to the universality class 
of the purely forward scattering model considered in [|r^. Indeed, in this special case the 
only field becoming massive is 9~ . In principal, it can't be ruled out that for some specific 
double chain models only one of two spin fields becomes gapful. However an investigation of 
the double chain Hubbard model [^ , demonstrates that the presence of the interchain one 
particle hopping is already sufficient to generate the antiferromagnetic spin exchange term 
with Jx ~ maxft u) ^hich will eventually make all spin modes gapful. We believe that it is 
a general feature of spin isotropic models of strongly correlated fermions on double chains. 
The author is indebted to Profs. T.M.Rice and F.D.M.Haldane for valuable discussions 
of these and related issues. This work was supported by the NSF Grant. 
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